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Abstract 

We generalize the predictions for attractions between over-all neutral sur- 
faces induced by charge fluctuations / correlations to non-uniform systems that 
include dielectric discontinuities, as is the case for mixed charged lipid mem- 
branes in an aqueous solution. We show that the induced interactions depend 
in a non-trivial way on the dielectric constants of membrane and water and 
show different scaling with distance depending on these properties. The gen- 
erality of the calculations also allows us to predict under which dielectric 

conditions the interaction will change sign and become repulsive. 
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I. INTRODUCTION 



In recent years there has been a growing interest in electrostatic systems that are dom- 
inated by ion fluctuations and ion distributions around larger charged objects. In some of 
these systems one finds attraction between like charged objects [0] and direct electrostatic 
contributions in systems that are over-all neutral 

In this paper we will generalize some theoretical results for systems of neutral surfaces 
(membranes) that nontheless interact electrostatically via ion fluctuations and correlations. 
These predictions are relevant to the experimental work done both on biological systems and 
on artificial systems where charges are introduced in order to improve membrane character- 
istics. Examples are the charged membranes in membrane-DNA complexes @] used for gene 
transvection and the formation of equilibrium bilayer vesicles from mixed charged lipids || . 

Recently it has been shown || that charge fluctuations can lead to attractions between 
over-all neutral surfaces. However, the system treated was the somewhat artificial case of 
uniform layers where the interacting surfaces separate regions of the same dielectric. In this 
paper we specifically focus on the role of the dielectric discontinuities in systems of lipid 
membranes in an aqueous solution and how they affect these interactions. In Sec. [II] we 
introduce a model system for the membrane which includes two surfaces charged with both 
positive and negative mobile ions (charged lipid heads at the bilayer surface) that are over-all 
neutral. The system is treated within the Debye-Hiickel model ||,[7| for a two-dimensional 
salt solution We calculate the interaction between these two surfaces resulting from 

the fluctuations and correlations of the mobile charges, and find that the resulting attraction 
depends in a non-trivial way on the dielectric discontinuity between lipid and water. 

II. INTERACTION BETWEEN TWO SALTY SURFACES 

In this section we calculate the effective interaction between two surfaces that contain 
mobile charges but are over-all neutral. This is a model system for mixed charged lipid 
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membranes |J or for membranes that are very highly charged to the extent that their 
counter-ions are restricted to a near-by layer that is thin enough to be considered as a two- 
dimensional surface. Pincus and Safran [Q] have calculated this interaction within the Debye- 
Hiickel approximation for a uniform system, i. e. a system with no dielectric discontinuities. 
We will follow their method, while introducing the dielectric contributions to this model. 



A. Model 

The Debye-Hiickel model is an expansion of the energy to second order in the charge 
density fluctuations [[UJ and includes both the electrostatic and entropic contributions due 
to these fluctuations: 



5H = Jdpdp' 



1 i=l,2 V ^0 / 

(2.1) 



The self energy of each of the surfaces separately is given by the first two terms while 
the third term is the interaction term between charges on the different surfaces. are the 
charge densities on the surfaces (the index % — 1, 2 denotes the surface number) while p is the 
in-plane coordinate and z is the coordinate perpendicular to the surface. The first term (5 
function) is the entropic contribution from the charge density fluctuations in both surfaces. 
In this expression we have assumed, for the sake of simplicity and without taking away 
from the generality of the treatment, that the charge fluctuations are only due to density 
fluctuations of one type of charge while the other sign does not fluctuate and therefore 
does not contribute to the free energy to this order. Thus the entropic contribution can 
be written in terms of the total charge density fluctuations on each surface , where <7o is 
the average charge density of each species (separately). The electrostatic contributions, 0, 
both between charges in the same surface (z = 0) and between charges on the opposite 
surfaces (z = d) are not trivial because of the dielectric discontinuities that are formed by 
these surfaces (fig 1). The discontinuities reflect the fields thus creating image charges in the 



region outside the membrane || . Because this system has two such discontinuities on either 
side of the membrane, each image charge is reflected over and over again so that we have 
an infinite number of charges over which to sum when calculating the potential. We require 
expressions both for the interaction between charge fluctuations in the same surface (they 
will also contribute to the inter-surface interaction via the reflections) and fluctuations on 
opposite surfaces. The interaction potential of two charges that are in the same surface is: 

0(p _ p", z = 0) = - ( — ^— + u 2 ^ 1 , 1 1 (2.2) 



While the interaction energy for two charges on the two different sides of the membrane 
is given by: 

4>(p- ft, z = d) = ^ f; u 2n ~ 2 . 1 (2.3) 

e n=i J\p-pf + ((2n-l)dy 



Here t wm are the dielectric constants of water and membrane lipid respectively, e 



e™±e™. ; u — e rn . <^w anc i ^ j g m embrane thickness. 



The sums in Eq. p.2| and |2.3| are easily performed if we use the identity / exp -92 Jo(qr)dq = 
-frTTp to transform them into simple geometric series. The resulting energy in momentum 
space has the form: 

1 



5h=y: 
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5aMV + IM?)r A (q) + 8<ri(q)6a 2 (q)B(q) 



(2.4) 



The coefficients are: A(q) = (± + ^f^ + ^f j^^^d) ) and B{q) 



. Here we have defined three different "Bjerrum lengths": <l>- 



2nl m exp(-qd) 
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At this point it is worth noting the differences between this expression and that which 
is found for the uniform case |2j] of no dielectric variations: The differences are expressed 
through the various effective Bjerrum lengths. In the uniform case there is only one such 
length scale which would be equal to <l> where e — e. In that case l m —<l>— Ib and 
51 = 0. Hence the differences enter not only in the way they change the interaction amplitude 
through l m and </>, but also by adding an additional interaction term that is d dependent, 



but which is also proportional to the dielectric difference, e m — e w , through 51, and affects 



the resulting interaction in a non-trivial way, as will be seen below. [ OJ] 

The Gibbs free energy for these fluctuations is now given by the logarithm of the partition 
function: 

^ = -log {| U q da q eM-^H/k B T)^ = log {.4(g) 2 - B(q) 2 } (2.5) 

The pressure between the two surfaces due to charge fluctuations as a function of mem- 
brane thickness is given by the negative derivative of the Gibbs free energy with respect to 
the thickness: 

A q (Xq + 1 + ^ exp- 2 9 d J - MjO exp- 2 <? d 

where we have introduced a Guoy-Chapman length scale: A = 27T <i>a • ^ n integral form we 
find the expression: 

TTfJ\ , rp 1 fi 2 -2x <S \d X + 1 + <Z> eXP ) ~ (<!>) / ~\ 

11(d) = kBT-^p / etaarexp za — ^ ^ — ^ ( 2 -7) 
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B. Results and Discussion 

The most convenient way to analyze the results of the previous section is by looking at 
the various limits of the integral, Eq. |2[7]. We have three dimensionless parameters that 
determine the behavior of this integral and thus the d and e w , e m dependence of the pressure. 
The first is the ratio between the two length scales in the problem: 

A _ 1 

d 2tt <l> ad' 

which parameterizes the strength of the charging in the membrane relative to the distance 
between the surfaces. The other two parameters are ratios of the dielectric constants and 
also their relative difference: 

51 2i(€ m CtoJCm i I'm 

and 
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The first of these two ratios reflects the effect of image charges on the fluctuation induced 
interactions, while the second measures the relative weakening or strengthening of the pri- 
mary interactions between fluctuations on the two sides due to the difference in dielectric 
response of the material between them. 

We have three different parameters with which we find two main regimes and one sub- 
regime. The first regime is reached when we take the limit ^ ■€ 1 (high ion density: the 
average distance between ions <C Vd </>): 



_, A k B T 51 ( l m \ \ k B T 
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The 1/d 3 behavior remains the same throughout this regime, although the sign of the pres- 
sure changes from being attractive for e w > e m (as is expected for a lower dielectric between 
the surfaces and is the case for a biomembrane) and even e m slightly bigger than e w , becom- 
ing repulsive only when the internal dielectric, e m is at least twice as big as the external one, 
e w . In this limit the effect of the variation in the dielectric between the surfaces is just on the 
size (and eventually the sign) of the pressure, but the dependence on distance is unaltered 
from the uniform case which was described in [Q] as a fluctuation effect and compared with 
the Van der Waals attraction also because of its linear dependence on temperature. 

The next main regime is the opposite one when 4 ^> 1. Here we distinguish between 
two regimes: The first is that when the dielectric contrast is not very big (compared with 
(<?>) anc ^ * n ^ ms case ^ e behavior is, as expected, similar to that found for the uniform 
case in this limit [01: 



lm V k R T fe m \ 2 cr 2 e 4 



In this case the pressure is inversely proportional to the temperature (through the A depen- 
dence) and is argued to be a correlation, rather than a fluctuation, effect [Q]. The dielectric 
effects enter in the coefficient (<?>) an d reduce the interaction as the internal dielectric 
(lipid) becomes smaller then the external one (water) and the dielectric contrast increases. 
However as this contrast increases another effect becomes important: the effect of the image 
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charges which dominate when |^^| is not small compared with y^f^j j, and we find: 

T-r/n SI k B T (e w - e m )e m ae 2 

n(rf)0C ^^A " P (2 ' 10) 

Here once again we find that the interaction will change sign when the internal and external 
dielectrics reverse roles. However, the dominant effect is that the power law changes from 
d^ 1 to d~ 2 , and therefore for smaller d this effect becomes more important than the previous 
result, Eq. |2.9| . Note that in this regime the pressure is independent of T and is therefore 
neither pure fluctuation nor correlation effect. Moreover, it is linearly dependent on the 
surface charge density, a, (and not quadratically) indicating that the correlations lead to an 
average charge distribution which is temperature independent and the result is an interaction 
between each charge and its effective image charge which does not include, to first order, 
the rest of the mobile charges. 

Because the membrane thickness is typically of order 40A the limit of 4 3> 1 can only 
be achieved for very low charging of the membrane and in this limit it might not be strong 
enough to compete with the Van der Waals interaction; in any case the stronger power 
dependence on d might not be easy to detect. However, note that if we reverse the membrane 
and water roles, and we look instead at the inter-membrane interactions say in a stack, we 
find that this last result might be more important. Because the inter-membrane distances 
in stacks can be relatively small this limit is easily achieved even in moderately charged 
membranes (one of every 5-6 lipids is charged). Because the dielectric constants are now 
reversed, the ratios change but we remain in this last limit where the reflections dominate 
the interaction. Moreover, due to the reversal of the dielectrics, the interaction (between 
membranes across the water) is repulsive and therefore the interplay with the Van der Waals 
attraction becomes more interesting. It is especially meaningful in this case because unlike 
the lipid material in the membrane, in some experimental set-ups water can flow in and out 
of the stack and therefore the stack separation can be more effectively controlled by this 
interaction. 

In summary, we have shown that fluctuation induced interactions are strongly dependent 
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on the dielectric properties of the system not only quantitatively but also qualitatively. The 
lower dielectric constant of lipid will reduce the strength of the interaction between the two 
surfaces of the membrane but will also change the scaling with the membrane thickness. 
When looking at interactions in a stack the reverse happens: the interaction is enhanced by 
a factor of ~ 2 with respect to the uniform case and we might also be able to see the 
effects of dielectric reflections when looking at the inter-membrane interactions. 
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FIGURES 

FIG. 1. Schematic of model system of membrane (e = 2) in water (e = 80). The lipid heads 
are charged both with positive and negative charges but the membrane is overall neutral. The 
dashed lines indicate the virtual surfaces where the image charges show up. Because there are two 
dielectric discontinuities, there are infinitely many such surfaces at equal distances, d, apart. 
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